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'. Abstract 

^ I We classify AqcD/w-b power corrections to nonleptonic B M1M2 decays, where Mi_2 are 

' charmless non-isosinglet mesons. Using recent developments in soft-collinear effective theory, we 

Q ■ prove that the leading contributions to annihilation amplitudes of order as{mb)AQCD /nib are real. 

. The leading annihilation amplitudes depend on twist-2 and the twist-3 three parton distributions. 

\^ I A complex nonperturbative parameter from annihilation hrst appears at 0[a'^{\^Kmb)AQC]:)/mi,]- 

I "Chirally enhanced" contributions are also factorizable and real at lowest order. Thus, incalculable 



• strong phases are suppressed in annihilation amplitudes, unless the as{\/ Ami,) expansion breaks 

down. Modeling the distribution functions, we find that (11 it 9)% and (15 it 11)% of the absolute 
Oh. values of the measured B'^ — > K^tt^ and B^ K~ penguin amplitudes come from annihilation. 

r~| ! This is consistent with the expected size of power corrections. 
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I. INTRODUCTION 



Nonleptonic charmless B decays are important probes of the standard model. They 
are sensitive to the CP violating phase 7 (or a) via the interference of tree and penguin 
contributions, and to possible new physics that could modify the penguin amplitudes. They 
also provide a powerful laboratory to study strong interactions, the understanding of which 
is crucial if one is to claim sensitivity to new physics in these decays. 

The theory of nonleptonic B decays underwent important progress in the last few years. 
Factorization theorems for B — > MM' decays have been proven to all orders in a.^ at leading 
order in A/mf,, for decays when M is a light (charmless) meson and M' is either charmed 
or charmless |]T], ^ |^, ^, ^ . Here A ~ Aqcd ~ 500 MeV denotes a typical hadronic scale. An 
important difference between the various approaches to making predictions for the charmless 
B M1M2 decay rates |, |, |, 0, |, |, 0, H HI' M is how certain 0{A/mb) power 



suppressed corrections are treated. In particular, it was observed that so-called annihilation 
diagrams (as in Fig. H) give rise to divergent convolution integrals if one attempts calculating 
them using conventional factorization techniques 0. In the KLS (or pQCD) approach [0, 
these are rendered finite by kj_ dependences, which effectively cut off the endpoints of the 
meson distribution functions. KLS found large imaginary parts from the jet scale, y/mi,A, 
from propagators via lm[xml — k'j_ + ie]~^ = —n6{xml — kj_) |jl4|. They also found that 



for the physical value of the power suppression of these terms relative to the leading 



contributions was not very significant. In the BBNS (or QCDF) approach p, 0, |TT[, the 
divergent convolutions are interpreted as signs of infrared sensitive contributions, and are 
modeled by complex parameters, Xa = Jq dy/y = (1 + pAe*"^-*) ln(mB/A), with Pa < 1 and 
an unrestricted strong phase ipA- In Ref. annihilation diagrams were investigated in the 
soft-collinear effective theory (SCET) |TB[ and parameterized by a complex amplitude. When 
annihilation is considered in SU{3) flavor analyses a complex parameter is also used |T^. In 
the absence of a factorization theorem for annihilation contributions, a dimensional analysis 
based parameterization with A/m^ magnitude and unrestricted strong phases is a reasonable 
way of estimating the uncertainty. In order not to introduce model dependent correlations, 
a new parameter could be used for each independent channel. 

It was recently shown by Manohar and Stewart [jl8| that properly separating the physics 
at different momentum scales removes the divergences, giving well defined results for convo- 
lution integrals through a new type of factorization which separates modes in their invariant 
mass and rapidity. The analysis involves a minimal subtraction with the zero-bin method 
to avoid double counting rapidity regions, and with the regulation and subtraction of diver- 
gences for large p"*" and p~ momenta that behave like ultraviolet divergences. Additional 
subtractions would correspond to scheme dependent terms, so the minimal subtraction is 
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the usual and simplest choice. We refer to this as MS factorization.^ In this paper we 
classify annihilation contributions to B ^ M1M2 decays and demonstrate how this rapidity 
factorization works for the leading terms of order 0[as{'mh)h./mb\. These leading order anni- 
hilation contributions are real despite the presence of endpoint divergences. We also classify 
which terms can involve a nonperturbative complex hadronic parameter, and show that they 
first show up for annihilation at higher order in perturbation theory, 0[a1{\/ mi,A) A/mf,]. 

Our analysis demonstrates that while certain annihilation contributions are only sensitive 
to the hard short distance scale /i^ ~ (local annihilation), there exist other annihilation 
contributions that start at the same order in as and l/rrib and are sensitive to the inter- 
mediate scale /i^ ~ mf,A (hard-coUinear annihilation terms). The leading local annihilation 
terms involve fs and a modified type of twist-2 distribution functions, while the leading 
hard-collinear terms have twist-3 meson distributions. In this work we perform matching 
calculations for the two-body distributions that require rapidity factorization. The calcu- 
lation of the leading amplitude involving the three body functions is given in a separate 
publication |jT9|, however we review the numerical results here. 



An interesting set of power corrections are those proportional to /ip where /i,r = ^^/{'mu + 
rrid) and /ix = m\/{mu + tRs) ||2^. For kaons and pions /ip ~ 2 GeV, so corrections 



proportional to fip/mi, can be sizable, and were labeled "chirally enhanced" in Ref. 0, [T3]. In 
the chiral limit /ip oc A^, where A^ is the chiral symmetry breaking scale, so the enhancement 
is not parametric, and comes from the fact that A^ > Aqcd- In the BBNS approach these 
/ml annihilation power corrections are included along with the leading order terms, and 
when they multiply divergent convolutions they are described by complex parameters. Below 
we show that, much like the lowest order annihilation contributions, these terms are also 
real and factorizable. 

In section |I| we review the leading order factorization theorem, and classify power correc- 
tions to -B — ^ M1M2, with a focus on annihilation amplitudes. In section |ITl| a factorization 
theorem is derived for local annihilation amplitudes at order A/mf, for final states not in- 
volving isosinglets (given in Eq. (|23|)). These amplitudes start at 0{as{rni,)) and involve 
and a modified type of twist-2 meson distributions. The extension to chirally enhanced local 
annihilation terms is considered in section |^ In section we study annihilation amplitudes 
from time-ordered products, and classify complex contributions generated at the hard scale 
nib, the intermediate scale \/mbA, and the nonperturbative scale A. Our results give abso- 
lute predictions for the annihilation amplitudes in B ^ PP, PV, VV channels, given the 



meson distribution functions as inputs, which are studied in Section This section also 
discusses the implications of our results for models of annihilation used in the literature, 
and a numerical analysis of the annihilation amplitudes in B ^ Kit and B — » KK. Ap- 
pendix gives the derivation of a two-dimensional convolution formula with overlapping 



^ Over the objection of one of the authors. 
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zero-bin subtractions. 



II. ANNIHILATION CONTRIBUTIONS IN SCET 

We use M to denote a charmless pseudoscalar or vector meson (vr, p, . . .). The 
relevant scales in B —y M1M2 decays are my/, rrify, E ^ "^5/2, rric, the jet scale vEA, 
and the nonperturbative scale A. Here E is the energy of the light mesons, which is much 
greater than their masses, mMi,2 ~ ^- To simplify notation, we denote by nib hereafter 
the expansion in all hard scales, {mf,, E,mc}. The decays B — > M1M2 are mediated by the 
weak AB = 1 effective Hamiltonian, which has AS* = terms for b —>■ dqiq2 transitions and 
= 1 terms for h —>■ sqiq2- For AS = it reads 

^ 10,77,89 

= 4 E ^p^K*^ (^1^1 + ^^of + Yl c'.o,) , (1) 

p=u,c i=3 

where the operators are 





= {ub)v-A {du)v-A, 




— (W/3&a)v-A {daUp)y_A , 


01 


= {cb)v-A {dc)v~A, 


01 


= (C/3&a)v-A {daC/3)v-A , 


03 


= T.q'{db)v-A{q'q')v-A, 


04 


= J2q'idf3ba)v-A {q'aq'/3)v-A , 


05 


= ^g'{db)v-A{q'q')v+A, 




= Y.q'idf3ba)v-A {^o^q'p)v+A , 


Or 


= Eg' 2 {db)v-A{qq)v+A, 


Og 


= Eg' 2 (dpba,)v-A{,qaq(3)v+A, 


O9 


= Y.q' 2 {db)v~Aiqq)v-A, 


Oio 


= J2q' 2 idpba)v-A{qaqp)v-A, 


O77 


= -^mbd<j^'F^,{l+^^)b, 




= -£^mbda''^G%T%l+^,)b 



Here current-current operators, a and (3 are color indices, Os-q are penguin 

operators and Oj^iq are electroweak penguin operators, with a sum over q' = u, d, s, c, b 
flavors, and electric charges Cg'. Results for AS* = 1 transitions are obtained by replacing 
d — > s in Eqs. (|1]) and (0), and likewise in the equations below. The coefficients in Eq. (|I]) 
are known at NLL order |]2T| (we have Of ^ Of relative to |^). In the NDR scheme, taking 
as{mz) = 0.118 and rrib = 4.8 GeV, 

Oi_io(mb) = {1.080, -0.177, 0.011, -0.033, 0.010, -0.040, 

4.9xl0~\ 4.6xl0-\ -9.8xl0"^ 1.9x10"^}. (3) 

To deflne what we mean by annihilation amplitudes we use the contraction amplitudes 
Ai, A2, P3, P^^^ in the full electroweak theory from Ref. (which thus includes pen- 
guin annihilation). These amplitudes are scheme and scale independent and correspond 
to Feynman diagrams with a Wick contraction between the spectator flavor in the initial 
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state and a quark in the operators Oj. Using SCET these annihilation amphtudes can be 
proven to be suppressed by A/rrif, to all orders in |]5|- These contributions differ from 
emission-annihilation amplitudes, EAi and EA2, which involve at least one isosinglet meson. 
As demonstrated in Refs. |]TT], EA12 occur at leading order in the power expansion. We 
focus on isodoublet and isotriplet final states, so ignore the -EAi 2 amplitudes hereafter. 

To separate the mass scales occurring below nib we need to match Hw onto operators in 
SCET. The nonperturbative degrees of freedom are soft quarks and gluons for the i?-meson, 
n-collinear quarks and gluons for one light meson, and n-collinear fields for the other light 
meson, as defined in [0]. Expanding in A/m^ gives 



{M,M2\Hw\B) = + A,, + 41 + + 



Gp^BfMi /a/2 fs 
V2A0 



A^^^ + Acc + A^n}„ + ^' 



(1) 

rest 



(4) 



In the second line we switched to dimensionless amplitudes A by pulling out a prefactor 
with the correct A^/^m^^ scaling. Here Aq = 500 MeV represents a i?-meson scale that is 
(9(Aqcd)- Taking r] = A/nib we have the leading order amphtude A^'^^ = 0{7f), and the 
subleading amplitude A*^^^ = Amn + Arrest = ^(^^)) which we have split into the annihilation 
amplitude A^ann and the remainder A^^}^^. The amphtude A^c in Eq. (^), denotes contributions 
from long-distance charm effects in all amplitudes, while perturbative charm loops contribute 
in the amplitudes A^^'^ and A^^\'^ 

There are two formally large scales, m?, ^ ^/miA ^ A, which we will refer to as the hard 
scale ~ fT^bi and intermediate or hard-collinear scale /Xj ~ y/rrihA. These scales can be 



integrated out one-by-one [0 with effective theories SCETj and SCETn. Integrating out 
TUh requires matching the Oi onto a series of operators in SCETi, Q^^^ ~ A-' where the SCETi 
power counting parameter A = rj^^'^ = A/A/m^. To obtain contributions to B ^ M1M2, we 
require an odd number of ultrasoft (usoft) light quarks q^s, two or more n-collinear fields, 
and two or more n-collinear fields, where = = 0. 

We briefly review results from Refs. [^, |^ for the leading amplitude A^^'^ for B — > M1M2. 
Here we have weak operators Qf]_Qd ~ A^, Qi^Lgd ~ -^^ with no gns's, taken in time-ordered 
products with an usoft-colhnear quark Lagrangian, C^^^ ~ A-' for j = 1,2, which has one Qus- 
We denote other subleading Lagrangians by C^^\ and hst the O(A^) and O(A^) time-ordered 
products for A^^^ in Table |l[ Matching these time-ordered products onto SCETn gives the 
leading 0{r]^) operators.^ When combined with the r/"''/^ from the states this yields a matrix 



^ Acc has the c- fields in Of 2 and O3-10 replaced by nonrelativistic fields ||], and is suppressed by at least 
their relative velocity, v ~ 0.3 — 0.5. The possibility of large nonperturbative charm loop contributions 
was first discussed in Refs. |l^, and the size of these terms remains controversial |2^, |2^ . 

^ Recall that to derive the rj^, we note that = V^, and changing the scaling A ^ 77 for four coUinear quark 
fields in matching SCETi SCETn gives the extra 77^. The A'^ term gains an extra A from the change 
in scaling to a coUinear D±. 
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element of order r/^/^, in agreement with the prefactor in Eq. @. Examples of the weak 
operators in SCETi are 



Qld = [^n,u>iigl$n,u^^PLhv] [^.t^a^^L^n.t^g] , (5) 

where other Qf^^^^ have different flavor structures. The "quark" fields with subscripts n and 
n contain a collinear quark field and Wilson line with large momenta labels cUj, such as 



[^^)Wn6{uo-n-r^)\. (6) 



Here creates a n-coUinear quark, or annihilates an antiquark, Wn = W[n ■ y4„] is the 
standard SCET collinear Wilson line built from the n component of n-collinear gluons, 
n ■ V'^ is an operator that picks out the large n ■ p label momentum of the fields it acts 
g, and igB^i;;, = [l/VWl[in-D^^n,iD'^^^]Wn5{uj - V^)]. The K is an HQET 6-quark 



on 
field. 



The leading order factorization theorem from SCETi is 



V2 



[ dudzTu{u,z)Cj^\z)(t)^'{u)+ [ duTx{u)C 
Jo Jo 



+ {Mi ^ M2}. 

(7) 

Here Tu and contain contributions from the hard scales mi,, and 0^ is the nonpertur- 
bative twist-2 light-cone distribution function. The terms (^^^ and Cj^\z) contain contri- 
butions from both the intermediate scale fii ~ i/mfeA and the scale A, and are defined by 
SCETi matrix elements between B and M states. In particular their scaling is 

C^^iE), Cr{z,E) ~ (^T + ...], (8) 

explaining the as{^J-i) entry in the A^^^ rows of Table § The C^*^ functions occur in both 
semileptonic decays and nonleptonic decays {E ^ 771^/2). Integrating out the scale Vm^A 
to all orders in as by matching onto SCETn gives [El ITS 



Cr E) = ^^^^ jdxj dk^ J{z, X, E) <p'^{x) <P^{k+) , 
C^{E) = ^^j^ E / 1 dktdkt Jabix,, kl, E) cp^ix.) 0f (A;;) , (9) 

a,b 

where the 0^ and 0f 's are twist-2 and twist-3, two and three parton distributions and we 
pulled out fsfu for convenience. The jet functions J, Jab occur due to the time-ordered 
product structure in SCETj and contain contributions from the scale y/rribK. Using the 
result for Cj^^ at order as{^i) this result agrees with Ref. (where expressing (^'^^ in 
terms of the full theory form factor generates an additional (f^^ term). The result for C^*^ 
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Order in 


Time-ordered products 


Perturbative 


order 


Dependence 


Properties 




m buijli 


Annihilation 


(Jtner 


m bUijin 


yl(0) 






as(/^j) 


(t)f (t)f 4)"'^ 


Real 












Real 


yl(l) 










Complex 




Qr 


as{^J'h) 






Real 




'-S_q 








Real 




l^^q J ' l-^^q J ^ 


al{^i)/TT 


a2(^i)/7r 


^ji^t,2^ ^3 ) ■ ■ ■ 


Complex 




^(0)r^(l)n2^(2) ^(l)r^(l)n3 


a^(^j)/7r 


a^(^i)/7r 


Ojl^fi,]^ 2) '*'3 / • • • 


Complex 












Complex 






as{^JLh) asifii)/TT 






Complex 


^(2) 








Real 



TABLE I: All contributions to i? ^ M1M2 amplitudes at leading order (A^^^) and at order A/m^ 
(A^^^), besides Acc- In the first A^^^ line / + j + — ^- The terms with — are absent or 
higher order when matched onto SCETn- The dependence in SCETn column lists the known 
dependence on nonperturbative parameters. The properties column shows whether at least one of 
the nonperturbative parameters is complex. For A^'^\ suppressed by A^/m^, only the local chirally 
enhanced annihilation operator is shown. 



is from Ref. 0] and required the MS factorization with zero-bin subtractions. The set of 
contributing functions (indices a, b) is determined by the complete set of SCETn operators 
derived in Ref. [Q. The power counting in as{fii) for the SCETj functions and (f^ 
agree with that derived in pQCD p9[] . 

Next we classify the contributions to the power suppressed B M1M2 amplitudes A^^\ 
In SCETi we need to study operators and time-ordered products with scaling up to 0(A^°). 
These have one or three light usoft quark fields. The relevant terms are listed in Table |, 
where ~ X^+^ and our notation for the Lagrangians up to second order is taken from 
Ref. [^. All the listed terms have an odd number of soft light quark fields. A basis for 
the Qf'^ operators is constructed in section |T|, for the Qf^C^^^ terms in Ref. |T^, and for 
the Q^^ terms in section |rV|. A basis is not yet known for the remaining qI^^ operators, 
for Qf\ and for the C^^^^^ and Lagrangians, but they do not contribute at 0{as), and 
only general properties of these operators are required for our analysis. Dashes in Table | 
indicate terms that are absent to all orders in for reasons to be explained below. To 
determine the perturbative order listed in the table we count the number of hard as{^h) 
factors from the matching onto SCETj, and the number of intermediate scale as{fii) factors 
from matching onto SCETn. The dependence in SCETn column lists the nonperturbative 
quantities that appear in the factorization theorem for the leading order result described 
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above, and from the factorization theorems we will derive in sections |TT| and |^ below. The 
properties column lists whether the nonperturbative distribution functions are complex or 
real as described in detail in section |V|, and has implications for strong phase information 
in the power corrections. The results in Table | imply the following power counting (for 
amplitudes not involving Acc), 

Re ~ asii^i) , Im ~ , 

nib 

RefiJiy ~ — , Im[iiiL] ~ "^(/^i) — • (10) 

nib TTT'h 



To facilitate the discussion we divide the annihilation amplitudes into local annihilation 

.(1,2) 

Lann 



contributions, A^l'^^^ from the operators Qf'^^ that are insensitive to the jet scale, and into 



the remaining annihilation amplitudes, ^^"ann' which are from time-ordered products in 
SCETi. Thus, 

Ail) _ Ai) , Ai) 



In the literature [0, ^, |TT], |^ only local annihilation amplitudes have been studied, 
and their matrix elements were parameterized by complex amplitudes. In SCET, Q^^^ is a 
six-quark operator with one usoft quark, such as 

where other qI^^ operators have different flavor structures. To derive the power counting 
for this operator, start with Q^^^ ~ A^, then note that switching a collinear quark to an 
usoft quark costs A^, and adding a and from a hard gluon also costs A^. This yields 
Q['^^ ~ 0{as{fih)^^^)- In matching onto SCETn we simply replace Qf^^ —>■ O^^^^ ~ rj'^, 
with the operator having an identical form. SCETi operators Qf'^ that do not have the 
form in Eq. (|1^) exist, but they must be taken in time-ordered products with a subleading 
Lagrangian and so do not contribute to A^^\ For this reason all local operator contributions 
to A'^'^^ contribute in the annihilation terms and not in A^^J^^. Since the matching onto 
O-^^'' is local, it appears as in Fig. ^ with an asi/J^h), but with no jet function. Thus 
this contribution to A^iln is of order a.sifJ'h) /o:s{f^i) -^/rrib relative to A^^\ In section ^ 
we construct a complete basis of Q-^^ operators and show that their matrix elements are 
factorizable in SCET at any order in perturbation theory, and do not generate strong phases 
at 0(a<j(/^/i))- We prove a similar theorem for chirally enhanced terms in the set Qf'^ in 
section 

The annihilation amplitudes and other A/nib suppressed amplitudes also occur through 
time-ordered products. Two examples are shown by Figs. |I|b and ||c. A subset of these terms 
were considered in Ref. [^], including the diagram in Fig. |l|c, and the phenomenological 



impact of these power corrections was studied. So far no attempt has been made to work 
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FIG. 1: Three types of factorization contributions to annihilation amplitudes which are the same 
order in r/ = AqcD/'Tih. a) shows Q^f^ which has > 1 hard gluon and factorizes at the scale 
rrih- The rapidity parameter, C, = p~ /p'^ , controls the MS-factorization between soft momenta 
(i?), n-collinear momenta (M2), and n-collinear momenta (Mi), b) shows the time-ordered prod- 
uct which involves factorization at and y^rn^. c) shows the time-ordered product 

Q^^^ [>C^g^]'^, which factorizes at the scale -v/m^A and does not need a hard gluon. Graphs a) and b) 
are of order as{nh)j while c) is as{ni)'^. 

out the strong phase properties and perturbative orders in as of the time-ordered products, 
a task we take up here. A complete classification of time-ordered products for the leading 
power corrections to i? ^ M1M2 is listed in Table |. A subset of these terms contribute 
to the annihilation amplitudes. To see which, we note that terms with a Qf'^^ and only 
one C^^^ do not contribute to annihilation at either leading or next-to-leading order; the 
weak operator is not high enough order in A to contain an extra n-n pair, and there are 
not enough C^g's to produce the pair through a soft quark exchange. To rule out these 
terms it was important that we are not considering isosinglet final states, which receive 
emission annihilation contributions already at leading order. The term Qi^\jC^^gY does not 
contribute to annihilation because we find that all annihilation type contractions are further 
power suppressed when matched onto SCETn. 

Time-ordered products with either a Q^^-"^^ or with three C^g^s do contribute to annihi- 
lation. Examples of these two types are shown in Figs. ||b and |I|c. Compared to the local 
annihilation amplitude from Q['^\ only the time-ordered product Q^"^^ d^^^ contributes at the 
same order in as- To demonstrate this, note that for terms with three C^q's all graphs have 
at least two contracted hard-collinear gluons and so are Graphs with a Q^^'^-* 
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start with one asifJ^h), and will also have an additional as{fii) from a hard collinear gluon, 
unless it remains uncontracted in matching onto SCETn. The uncontracted gluon costs an 
additional A in the matching onto SCETn, so only the time-ordered product Q^'^'^C^^^ can 
have a leading, 0{as{mh)), contribution. Fig. gives an example of a diagram occurring 
from this time-ordered product. The resulting amplitude involves the three-parton distribu- 
tion, 03A/2- As shown in Ref. |T9[ it also involves the twist-2 distribution (^^, and its leading 
order convolution integrals converge. 

The time-ordered products with three C^g^s are suppressed by a'^{fii)/as{fih) relative to 
Qi^"* , and can be proven to involve a complex nonperturbative function, as labeled in Table | 
(an example is shown in Fig. |T|c). Thus, if perturbation theory converges rapidly at the 
scale /ij, then complex annihilation amplitudes are highly suppressed. If perturbation theory 
at /ij is poorly convergent then the time-ordered product contribution could be important 
numerically; comparable or even larger than the leading local annihilation amplitude from 
Local annihilation contributions are discussed in detail in sections |T| and while 
strong phase properties of the amplitudes and the time-ordered product contributions are 
taken up in section 0. 



III. LOCAL SIX-QUARK OPERATORS IN SCETn 

In this section we construct a complete basis of Of'^^ operators in SCETn (the Q-^-* 
terms in SCETj) and derive a factorization theorem for their contributions to B ^ M1M2. 
To find a complete basis we consider color, spin, and flavor structures that could appear 
when matching at any order in a^. Color is simple, the six-quark operator must have 

® r„ = l ® l ® l. Although operators with a in one or more bilinears are 
allowed at this order, with the factorization properties of the leading Lagrangians and 
{MnMn\0\Bs) = (0| . . . \Bs){Mn\ ■ ■ ■ |0)(M„| . . . |0), the terms with T^'s give vanishing ma- 
trix element between the color singlet hadronic states 

For spin we start by looking at chirality which is preserved by the matching at mh. Since 
there is no jet function, the soft spectator quark that interpolates for the S-meson must come 
from the original operator in H]^, and we Fierz this ip field next to the 6-quark field. To be 
definite, we take the other ijj field from Hw to go in the n direction (in the SCET Hamiltonian 
we sum over n ^ n). This implies that the pair-produced quark is in the n direction. For 
Oi-4,9,io the allowed chiral structures induced in SCET by matching are {LH){LL){LL) and 
{LH){LR){RL) where L and R correspond to the handedness for the light quarks in the 
bilinears in the order shown in Eq. (|T^ . We cannot assign a handedness to the heavy quark 
denoted here by H. For O5-8 we can have (LH) (RL) (LR) , (LH) (RR) (RR) , (RH) (LL) (LR) , 
and {RH){LR){RR). A complete basis of Dirac structures for the individual bilinears is: 

r, = 7", Tn = {i,^ll}, r„ = {^,^7^^}. (13) 
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Structures with 75 are not needed because we have specified the handedness. Here fl7^ 
and i/i'-)^^ connect left and right-handed quarks, while t/l and i/i connect quarks of the same 
handedness. From the basis in Eq. (|T3|) we must construct an overall scalar using the 
tensors f^, n^, n'^, g'^, and e'^ = e'^'^"^?7,Q,?T,^/2. We take e""*^^^ = 1, and work in a frame 
where f ^ = and n ■ v = n ■ v = 1, which makes the set {n, n, f } redundant. For reasons 
that will become apparent we pick and (n^ — ra^) as our basis in this section. The definite 
handedness allows us to turn any contraction involving ie^^ into a contraction with , for 
example le^^il^i^^i^ = ^f^i^ll^^^ = ^nillit The {LH){LR){RL) and {LH){RL){LR) 
structures can be ruled out since 

^7^^Pr ® ^7^ Pi = ^7^^Pi ® ^7^P^^ = . (14) 

Noting that = this leaves four allowed spin structures 

r, ® Ffi o r„ = {1 ® ^ ® ^, {^-^l) ®i®ifi, 7^ ® iii ® i] ■ (is) 

The last two structures have qsl'xpv and vanish identically for P-meson decays (they would 
contribute for P*'s). Furthermore, the local annihilation operators are not sensitive to the 
momentum of the soft spectator quark. Thus in taking the matrix element we can use 

(0|g,75/i.|P) = -imsfB , (0|g,75(^-?^)/i.|P) = . (16) 

Here is the decay constant in the heavy quark limit. The fact that we can match onto 
a basis of local SCET operators of the form in Eq. (0) demonstrates to all orders in as 
that the local annihilation contributions are proportional to ]b- Using Eq. ([I6D the second 
Dirac structure in Eq. (^) is eliminated, so we do not list operators with iifi—iji) in the soft 
bilinears below. 

Next we consider the allowed fiavor structures. From the operators Oi_2 "we have 
{ub){dq){qu), {db){uq){qu), from Oi-ej^.sg we have {db){q'q){qq'), {q'b){dq){qq'), and O7-10 
give a combination of these. Here the qq are the pair produced n and n pair, while the q'q' 
appeared in the weak operators. Thus a basis for P-decay operators is 

2 



m 



b 



2 

O^d^^ = —J2q [UsPRbv] [dn,u,2iPL Qn^u^z] [gn.i^i^^LMn.o;*] 



ml 



2 

m^ 
2 

Oid^ = — Eg,g' [Q'sPnbv] [dn,c.2iPLqn,c.,] [gn,<^i#^Lg^,^J 

2 

^M^^ = —J2q [dsPRbv] [Un,uJ2^PRqn,LU3] [qn,uJi'^PR Un,uJ4] , 



mr 



2 

C'ed^^ = — 3 Eg [UsPRbv] [dn,cJ2^PRqn,uJ3] [^n.o;! Mn,<^4] 
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2 

Here we integrated out c and b quarks in the sum over flavors, so the remaining sums are over 
q = u,d,s and q' = u, d, s. For the AS" = effective Hamiltonian with Wilson coefficients 
al'^\uj) we use the notation 

n.n 1=1—8 

To pull the CKM structures out of the SCET Wilson coefficients we write 

, f Xu'' aiu{ujj) + X^f^ aiciuj) [i = 1, 2, 3, 4] , 

[{X^u^ + Xr)ai{ujj) [z = 5,6,7,8], 

where X^^ = VpbV*^. Identical definitions for a| are made by replacing xlf^ —>■ •* and 
A^*^) ^ Xi'l For 2 = 5, 6, 7, 8 only penguin operators contribute. 

Next we take the B M1M2 matrix element of Hw- The factorization properties of 
SCET yield 

+ {Mi^M2}, (20) 

with similar results for the other O-]'^'' terms. Here the {Mi ^ M2} indicates terms where 
the flavor quantum numbers of the M2 state match those of the n-collinear operator. The 
matrix elements in Eq. (^OD are zero for transversely polarized vector mesons in agreement 
with the helicity counting in Ref. Equation (1^) can be evaluated using Eq. (|16D and 



{Pni{pWJl^PL,Rqn,2\0) = Cpff> 6nm S{fl- p -U + Uj') (ppijj) , 

{VnAP,eMliPL,Rq[{'M = '-^^^^^Cyff,6^^J{n.p-u+J)<Py^^{y). (21) 

Here /, /' are flavor indices, 4>p{y) and (py^^ iv) are the twist-2 light-cone distribution functions 
for pseudoscalars and vectors, y = uj/n ■ p = uj/rrib, and cpff, cyff are Clebsch-Gordan 
coefficients. For the M2 mesons, and Vn^-, we have the same equation with n n, 
and y ^ X. Since the Pl,r only induce ± signs in the pseudoscalar matrix element, it is 
convenient to define 

= + , ^2 = 02 + /tag , aj = a'^^ + Ka^ , 5^ = 04 + /tag , (22) 
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MiM2 




■K~TT~^, 7r~/9+, p~7r+, 

TT^TT^, p^TT^ T^ P^, P\\ P\ 
_^(*)0^(*)0 


-af (x, y) - dj{y, x) - of (x, y) - dj{y, x) 
[diix, y) + di{x, y) - x) - dfiy, x)] 

rij, , J, . 1/7/ r ~\ 

[iaf (x, y) + (x, y) + ^al{x, y)\ + [x ^ yj 
-of (x, y) - a^(x, y) - d'^{y, x) 
a;^(x, y) + a^(y, x) + a^(x, y) 
d^{x,y) + df{x,y) 


vr ' , p ' 


d2{x,y) + d%{x,y) 
-^[m^^y) +a|(x,y)] 
-^al{x,y) 
-a|(x,y) 



TABLE II: Hard functions for and B decays for the annihilation amphtude ^4^^^^^ in Eq. (| 
For each pair of mesons in the table, the first is Mi and the second M2. 



MIM2 


i7(x,y) 




-a|(y,x) 
■;75a4(y>2;) 


vr^vr^, vr^yO^, p^n^, p^ p^ 
i^(*)o_^(*)o 


-af(x,y) - a|(x,y) - a|(y,x) 
[iaf (x, y) + a|(x, y)] + [x ^ y] 
-of (x, y) - d|(y, x) - a|(x, y) - d|(y, x) 
y) + d|(y, x) + d|(y, x) 



TABLE III: Hard functions for Bg decays for the annihilation amplitude in Eq. (p^). 

—1 for PV channels. 



with similar definitions for a|. Here k = +1 for PP, V^V^, and k 
Using these results, the 0{K/mb) local annihilation amplitudes are 



M1M2 



V2 



dxdyH{x,y)(t)^\y)(t)'^\x) . 



LM21 



(23) 



Here H{x, y) are perturbatively calculable hard coefficients determined by the SCET Wilson 
coefficients di{ujj). Results for different final states are listed in Table |I| for P° and B~ 



decays, and in Table |T| for P^ decays. Our derivation of the local annihilation amplitude 
in Eq. (|2^) is valid to all orders in a^, and provides a proof of factorization for this term. 

Matching at tree level, involves computing the 0{as{mh)) graphs in Fig. ^and comparing 
them with matrix elements of the SCET operators Qf"^- Doing so we find that the Wilson 
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d) 



X 



X 



FIG. 2: Tree level annihilation graphs for B — > M1M2 decays. Here soft, n, n denote quarks that 
are soft, n-coUinear, and n-collinear respectively. 



coefficients ai{x,y) are 



O'Au 

a^ 



Ar2 

c 
c 

Cir7ra,(/ifc) 

c 
c 

CpT^asilJ'h) 

c 



F{x,y 
F{x,y 
F{x,y 
F{x,y 
F{y,x 
F{y,x 



c 



10 ; 



C2 + 2 ^9 



C4 ~ 2 



— — Co 

J 2 y 



Ar2 

c 

m 



F{x,y) (^Cio), 



( 2^9 



F{x,y){Cs--C, 



2 



1 



Og = . 
as = . 



(24) 



where fih ~ "^b? 3; = 1 — x, ?/ = ! — y, with quark momentum fractions x and y as defined 
in Eq. ( |2l] ) and shown in Fig. 0. The function F is 

rf(/i_) <5'(s) 



F(a;,?/) 



+ 



y 



(25) 



?/(x?/-l)_ 

where the 0-notation and term involving the Wilson coefficient d{fi-) are discussed below. 
The function F{y,x) will involve (i(/i+). Note that the coefficients a3u,3c,iu,4:c,7,8 are polluted 
in the sense of Ref. 0, meaning that 0{al) matching results proportional to the large 
coefficients Ci,2 could compete numerically. The others are not polluted: ai„,2u involve Ci,2 
at 0{as), while ai^ 2c,5,6 only get contributions from electroweak penguins. Our results for the 
diagrams in Fig. |^ agree with Refs. 0, |lOl . This includes the appearance of the combinations 
of momentum fractions in the functions F{x, y) and F{y, x), up to 0-distribution and d-term. 
For later convenience we define moment parameters which convolute the hard coefficients 
with the meson distributions 



P^^^''= / dxdy[aUx,y)+na,+,ix,y)]<P^''iy)^^''ix), 
Jo 

pmM,^ I c?a;rfy[a,,(x,i/)+/€a,+4(x,y)]0*''Ky)0*'Ha;)- 



(26) 
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In Eq. ( pSD the subscript denotes the fact that singular terms in convolution integrals 
are finite in SCET due to the MS-factorization which involves convolution integrals such as 

f dxr dx[. 6{l — x — x') '^^'^^ (27) 

1') 



where x^'-* and xf correspond to label and residual momenta [jl8|. Implementing x 7^ and 



x' 7^ in the MS-factorization scheme requires zero-bin subtractions and divergences in the 
rapidity must also be regulated. The 5-function sets x' = 1 — x, so x' 7^ enforces x 7^ 1. 
With the usual assumption that (pu^x) vanishes at its endpoints with a power-like fall-off 
slower than quadratic, only integrals over 1/x^ in F(x, y) and in F{y, x) require special 
care, 



= ^ . (,-^>" = j\ ^ . (28) 

'0 \^ )0 Jo [y )0 



The resulting moments {x'"^)^'^ and {y^"^)^^ should be considered hadronic parameters, for 
which we use the minimal subtraction scheme. Their value depends on /i and fi± and are 
scheme dependent beyond the usual MS scheme for 0m- This can be viewed as a modification 
of the distribution function, (pM^x, /i) (t)M{x, fi, fi-), where the x~^ moment of (pM^x, fi, /i_) 
converges. In order to derive a result that makes it easy to find a model for these moments 



we follow Ref. and assume there is no interference between the rapidity renormalization 



and invariant mass renormalization, which gives 

n X^ \ U- / 



X V /i 



Here 0^(1) is generated by a zero-bin subtraction which avoids double counting the region 
where x ^ 0. When x ^ the corresponding outgoing quark becomes soft, and this 
contribution is taken into account by a time-ordered product term in Table |. To obtain the 
renormalized (x~^)^ result in Eq. ( p9D requires 1/euv counterterms which correspond to 
operators with the ra-collinear bilinears in Eq. (p!7D, [Un,ui2i>'l5'in,ui-^ etc., which can be written 

m 



as 



Oct = -^i^nWhMW^^nU 



(30) 



The matrix element of these terms is taken prior to performing the partial derivative and 
the limit uj-^ — > 0, and gives 0'^^(l,/i). These terms do not have a tus 7^ restriction, and 
consistency of the renormalization procedure used to obtain Eq. (p9D demands that the 
fields here are n-collinear. An analogous set of terms are required for 0'jv/(O, /i). These terms 
are real at any scale, which follows from the requirements discussed in section for an 
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SCETii operator to be able to generate a physical strong phase. The dependences on /i-t in 
Eq. (|29| ) are canceled by the leading dependences on these scales, d{jj,J) = ln(p^^/yU_) + k 
and d{fi+) = ln(p^//i+) + k, which appeared in Eq. (|25|). Here k can be fixed by a matching 
computation. The d{fi±) correspond to the renormalized coefficients of the Oct, and must 
be included for consistency at this order In the rough numerical analysis we do later 
on, we will treat the contributions from these coefficients as part of the uncertainty. 

Note that in deriving the result in Eq. ( P^D we have dropped ie factors from the propa- 
gators. If these terms were kept, the second term in F{x,y) would be 

(31) 



{y + ie) {xy — 1 + ie) 



The ie's yield imaginary contributions with 6{y) and S{xy — 1). They contribute for y = 
or for X = y = 1, so these contributions occur in zero-bins, which are excluded from the 
convolution integrals in the factorization theorem we have derived with SCET. The zero-bins 
correspond to degrees of freedom that are soft, and including these regions would induce 
a double counting, so the correct factorization theorem in QCD does not include them. 
Factors analogous to x 7^ and x' ^ in Eq. (|27|) ensure that there is no contribution 
to the integral from any zero-bin momentum, and we find that the 5-function terms give 
zero. This remains true for more singular distributions yielding S^"'\x), and so also applies 
to the first term in F{x,y). Thus it is correct to drop the ie factors from the start. This 
should be compared with the approach in KLS where the ie factors generate a strong phase 
from the tree level diagrams from a kj_ dependent ^-function. In our derivation any such k'j_ 
imaginary terms could only occur at higher orders in A/mf,. 

Thus at order OsifJ^h) the lowest order annihilation factorization theorem is determined 
by the convolutions 



dxdyF{x,y)^^^'{y)4>^^%x) (32) 

= {x-'f\y-'f' - (b(xy- 1)]-^)^^^^ + rf(/._)0;,,(i)(r^>"'^ , 

1 

dxdyF{y,x)(l)^^^'{y)(l)^^^{x) 
= {y-'f\x-Y^ - {[xixy - l)]-Y'''' - dif,^)<t>'MM^-'f' ■ 



Here we use Eq. (pQ]), and 



{y-T= fdy^^^, {f{x,y)r^''^= fdxfdyf{x,y)ct^'Hy.l^)<P''Hx.l^)- (33) 

'0 y Jo Jo 



These results do not have a complex phase because the right-hand side of Eq. (Pz]) is real. 

We have shown that the convolution formula in Eq. ( p3D for the local contributions O-^^^ 
yields a well-defined annihilation amplitude. At order asinif,) the result is real, so is 
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real up to perturbative corrections. Order a'^{mb) corrections to the will produce pertur- 
bative strong phases in A^^j„„. Further discussion on strong phases is given in section |V], 
while phenomenological implications are taken up in section |VI[ 

IV. CHIRALLY ENHANCED LOCAL ANNIHILATION CONTRIBUTIONS 

At order as{Hh)f^M-^/^b there are contributions from chirally enhanced operators that 



could compete with the as(/i/i)A/mfc terms [|T0|- In SCET we define these contributions 
as the set of SCETn operators analogous to Of'^^ but with an extra between collinear 
quarks fields. We start by constructing a complete basis for local operators at this order 
with a V^, calling them Ol^^\ These operators have the same color and fiavor structures 
as Eq. (p!?!). The chiral structures induced from the operators Oi_io and the initial basis 
of Dirac structures shown in Eq. (|13|) are also the same, and allow us to eliminate many 
possibihties. 

The complete set of Dirac structures from matching the operators Oi_4_9,io include 

7^ ® ^ll ® 07^^^ , ll ® H ® ^^aK } , (34) 

plus the analogous set Ts^TnV^^Tn- Our basis does not include operators with V]_, because 
the mesons Mj have zero _L-momenta, so we can integrate these terms by parts to put them 
in the form in Eq. (p^ . The third term in Eq. (|3^) has chiral structure {LH){LR){RL) and 
vanishes by Eq. (0). The terms in Eq. ( p4D all have (jsl^by, and so do not contribute for 
5-decays. The same holds if we replace Vf_ by igB^. Thus, at any order in perturbation 
theory the only (9(?7^) local operator contributions from Oi_4^9^io are those with a in the 
soft bilinear. 

For O5-S we have the structures in Eq. (p^, and when the q' fiavor is a soft quark with 
Pl ® Pr Dirac structure from Oj we also have 

r, ® TnV^ = {1 ^ ® ^f^, 1 ® ^7^ ® } , 

Fs ® TnV^ ® r„ = {1 ® ® 07^, 1 ® i/if^ ® } , (35) 



plus operators with 1 replaced by ^ — which vanish due to Eq. ([iq) . The operators in 



Eq. (|35D contribute to 5-decays. In particular, they yield both transverse and longitudinal 
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polarization in B VV. A complete basis for the local 0{ri^) operators with one IS 



O 



(2L) 
Id 



^2d 



ml 



ml 



^Pr 



Csf^ = £4 WsPlK] [dn,u.2^1^PRqn,u.3] [qn,u;,^PRV^q'^^^^] 



o 



(2L) 
Ad 



o 



(2L) 
5d-8d 



ml 



ml 



Id-id 



(36) 



with sums over g, q' = u, d, s. Note that the flavor structure of these operators is identical to 
O^lP . For the the electroweak penguin operators Or^g an additional four operators OfJ^^^ 

(2L) 

are needed, which have the same spin-flavor structures as 0\^_^^, but with an e,/ charge 
factor, '^qqi 3eq'/2. Again we caution that we have not considered the complete set of local 
A^/m^ operators, since our basis does not include three-body terms with an igB'^, nor terms 
with an extra Dg soft covariant derivative. We have also not considered C(/XMi/^M2A/m^) 
terms. All these terms are real, and it would be interesting to calculate them in the future. 



(2L) 

The weak Hamiltonian with Wilson coefficients for the operators O-^^ is 



w 



^ (A(f) + Ai'^)) J2 1 [dooidoo2duj,duj4 J2 

* n,n 1=1— i 



(37) 



Since only the penguin operators Os-g contribute, we pulled out the common CKM factor. 

(2L) 

Matching at tree level onto the operators O^^ by keeping terms linear in the ±-momenta 
in Fig. 0, we find 



a^{x,y) 
a}{x,y) 



N'c 

4Cj7^7rag(/ifc) 
iCFTrasifih) 



{Ce + ^)Fr{x,y) + ^F,{x,y) 



a. 



c 

(c, + |)F3(..,)-|F.te.. 



05-8(3;, y) = a'^-iix, y) with C5 C7, Cg ^ Cg , 



(38) 
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where x and y are defined in Fig. 

l + x 



and 



Fi{x,y) 
F^{x,y) 



y'^ y 



+ d^{^^.)5'{x) 



y^y 



+ d2ifi+)6'{y) 



l + x 



x^ 



+ dsifi^)6'ix)6'iy) 



y2 ^2 



xy)xy'^ 

+ (i4(yU_)5'(x) 



r 



+ d^{ii+)5\y) 



x^ 



+ d,{ii±)5\x)5\y). (39) 



Here c/i_6 play the same role as d in Eq. (^). The coefficients polluted in the 

sense of Ref. meaning that 0{a1) matching results proportional to the large coefficients 
Ci_2 could compete numerically. This makes the computation of these 0{a1) corrections 
important. 



i2L') 

For decays involving a pseudoscalar in the final state, the operators and O 



){2L) 
'2d 



generate so-called "chirally enhanced" terms, proportional to /iM- Time-ordered products 
of SCETi operators also generate [im terms, but only at 0{a1). It is not clear that the 
chirally enhanced terms are larger numerically than other power corrections. In particular 
three-body distributions from operators with C,n{w^±)^^n are parametrically (and some- 
times numerically as well) of similar importance |3^. The distributions are related by p3[ 



fpfJ'P 



fpfJ'P 



kP'i 



x) + 



{2x - 1) 
x{l — x] 

1 



[X] - 



Qxil 



X) 



-6/: 



3P 



G 



X] 



G 



+ 



3P 



X] 



X 



1 — X 
1 — X 



(40) 



where Gp^(x) and G'p^(x) are integrals over the three-parton distribution, 03p. These re- 
lations allow certain chirally enhanced terms with iipfp to be traded for non-chirally en- 
hanced terms with /sp. Thus it is clear that the chirally enhanced terms dominate over the 
three-body operators only in the special case when the linear combinations in the square 
brackets on the left-hand side of Eq. ( ^D] ) are numerically suppressed. Solving with these 
linear combinations set to zero determines the two-body distributions 0^ and 0^ in the 
Wandzura-Wilczek (WW) approximation Thus in order to uniquely specify the /xp 



dependent terms, the WW approximation was needed in Ref. [T^ . 

In contrast, in SCET we are not forced to assume a numerical dominance of the /ip 
terms to uniquely identify them. We can instead define local chirally enhanced annihilation 



terms to be the matrix elements of the operators Of^^ and Og^^'' for final states with a 
pseudoscalar. With a minimal basis of operators, the matrix elements of these terms are 
unique. The remaining terms involve other operators, and we postpone discussing them to 
future work. We proceed to work out the factorization formula for 0[^^^ and 0^2d'^ with 
steps analogous to Eqs. (|20| ) through (pSj). To take the matrix element we need Eq. ( pT]) 
and the result 

-,U) 



,(2L) 



PnAvMJlif^PRq^C |0) = ^Pff ^nrn S{n-P -co + co') /p/ip0^,(y) . (41) 
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Here cp/// are Clebsch-Gordan factors, y = u/n ■ p, and we have not written the u' depen- 
dence in the distribution due to the (5-function. The distribution 4>pp{y) is related to more 
standard twist-3 two-parton and three-parton distributions by [|1^, Q 

2/3P f(¥ 

y' 



Note that in fip^pp, 



3y 



(y) 



hp{y-y',y) 



(42) 



'-^ the 03P term does not have the chiral enhancement factor /ip. There will 
be additional terms proportional to 03P generated by three-body operators. We choose the 



For 



(j)pp and (f)3P basis of twist-three distributions, keeping in mind the relations in Eq. 
decays involving one or more pseudoscalars in the final state we find the chirally enhanced 
local annihilation amplitudes 

,(2) _ GpfsfMifMi 



A 



Lann 



nib 



(Ai^) + A(^)) / dxdyLM,H^i{x,y) 0jJ(l/)0''H^ 



+ ^^M,H^2{x,y)<P^'Ky)<Pl'ix 



(43) 



where fip = fix* = and using isospin /x^ = m^/imu + m^), = ^k/{^s + n^^u) = 
m\/{'ms + md). Terms with 03P or terms of the same order with a D'^ in their soft matrix 
elements have not been included in our ^^.ann' though they also give local annihilation con- 
tributions to A^'^'^ . Furthermore, we focused on the pseudoscalar matrix element in Eq. (^Tj) 
to derive the contribution in Eq. (|i3|). The 0\^2d operators in Eq. ( p6D will contribute ad- 
ditional terms for decays to longitudinal vector mesons involving distributions h\^^'' and h^^^ 
(our notation for these distributions follows Ref. |0). The operators Of^l^ will produce 
decays to two transverse vectors with distributions from among F, V, A. It would be 
straightforward to work out a factorization theorem from the operators in terms of 

these distributions, though we will not do so here. 

Results for the hard coefficients H^i and H^2 in terms of the Wilson coefficients are 
given in Table |I^ for and decays and in Table for Bs decays. Note that there are 
no chirally enhanced annihilation contributions for the Bg or Bg —>■ pn channels, so 

Bg decays could potentially be used to separate annihilation contributions from A^^l^^ and 
^Lann- "^^^ later convenience we define moment parameters 



M1M2 



/0M1M2 



- dxdyal,{x,y)ct>f;{y)ct>''-{x) 



dxdyal,{x,y)<l)''^{y)<j)^l%x) 



(44) 



Neglecting 03P in the WW approximation yields (ppp{y) = 6?/(l — y)- At order as{ph) our 
results for P^i and P^2, taken with the WW approximation, agree with the convolutions 



derived in this limit in Refs. [Il0| , |TT|. Ignoring the 0-distributions we would find that these 
convolution integrals diverge. The zero-bin avoided double counting in our convolutions, 
and yields a finite and real result for the chirally enhanced annihilation amplitude. 
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MiM2 

TT^ir^, p^TT^ TT^ p^ 
TT+Vr^, TT^p^, p^TT^ 



7r-K(*)o, p-K^ 
■K+ K^*^~ , p+ K- 



Hxi{x,y) 

-^a^{x,y) --^a^{x,y) 
-^a^{x,y) + ^a^{x,y) 
-a^{x,y) + ^a^{x,y) 
a^{x,y) - iaf(x,y) 

ai{x,y) - ^a^{x,y) 
a^{x,y) +a^{x,y) 



ii{x,y) +a^{x,y) 
7^ai(.x,y) - -^a^{x,y) 



V2 

1 nX 



V2' 



i{x,y) 



■a^{x,y) 



-a^{x,y) + ^a^{x,y) 



Hx2{x,y) 

a2{x,y) - ^a^ix,y) 
-a^{x,y) + ^a^{x,y) 

-a^{x,y) + ^a^{x,y) 
-a^{x,y) - a^{x,y) 



-a^{x,y) - a^{x,y) 
-l=a^{x,y) + -i=a^(x,y) 

02 (3;, y) - ^(4{x,y) 



TABLE IV: Hard functions for the annihilation amphtude in Eq. ( ^3[ ) for and decays. 

The result for B^ tt^tt^ is obtained by adding the results using the entries from the first two 
rows, and so vanishes in the isospin limit. 



M1M2 


Hxi{x,y) 


Hx2{x,y) 


K+ir-, K*+TT-, K+p~ 


-ai{x,y) + ^a^{x,y) 
-j^a^{x,y) - j^a^{x,y) 


a^{x,y) - ^a^{x,y) 
-^a^{x,y) + ^a^{x,y) 


K+K-, K*+K-, K+K*^ 


-ai{x,y) + ^a^{x,y) 
ai{x,y) - ^a^(.x,y) 


alix^y) - \al{x,y) 
-a^{x,y) + \al{x,y) 



TABLE V: Hard functions for the annihilation amplitude A^^^^^ in Eq. ( ^3|) for Bg decays 



Let's see how the convolutions work out at order as(/i/i) following Ref. |T8|. We need two 
standard convolutions involving zero-bin subtractions, 
"1 

dx dy 



l + x 



dx dy 



y"^ yx"^ 

i + y 

y'^ xxP' 



Here we model the y~^, y~^ moments as in Eq. ( pQ]) and Eq. (|33D, and for the remaining 
convolution we again assume there is no interference between the rapidity renormalization 
and invariant mass renormalization to find 



/ _2 — 

\y y / pp 



dy 



+ 



PP 



(0,/i) In 



(46) 
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The /i-t dependence is canceled by tree level logarithmic dependence in the coefficients, 
rfi,4(/i_) = \n{pjjfi_), rf2,5(/i+) = ln(PM/^+), c?3,6(/i±) = ln(PA//^-)ln(p^//i+). The kernels 
in Eq. (^8[) also involve two more complicated convolutions that are derived in Appendix 0, 
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As promised, the minimal subtraction scheme yields a well defined result for A\l^^^. The 
scheme dependence cancels order by order in between the matrix element and perturbative 
corrections to the kernels obtained by matching. In any scheme the result at order asifJ^h) 
is real. 



V. GENERATING STRONG PHASES 



In this section we derive results for the order at which strong phases occur in the power 
suppressed amplitudes A^^\ It is convenient to classify complex contributions to the B — »• 
M1M2 amplitudes according to the distance scale at which they are generated. We use 
the terminology hard, jet, and nonperturbative to refer to imaginary contributions from 
the scales m^, y/rubA, and A^ respectively. We will not attempt to classify strong phases 
generated by charm loops, since a complete understanding of factorization for these terms 
order by order in a power counting expansion is not yet available. 

For a matrix element to have a physical complex phase it must contain information about 
both final state mesons. Generically, terms in the factorized power expansion of B ^ M1M2 
amplitudes involve only vacuum to meson matrix elements, so strong phase information can 
be contained in the Wilson coefficients or the factorized operators, but not in the states. This 
provides tight constraints on the source of strong phases. Nonperturbative strong phases will 
occur if matrix elements of these factorized operators give complex distribution functions. 
A sufficient condition to generate a nonperturbative phase, is to have a factorized operator 
that is sensitive to the directions of two or more final state mesons , information that can 
be carried by Wilson lines. Physically, this is a manifestation of soft rescattering of final 
states. In processes like ours where soft-coUinear and collinear(n)-collinear(n) factorization 
are relevant, and there is only one hadron in any given light cone direction, this criterion 
implies that all strong phases reside in the soft matrix elements, where the directional 
information from collinear hadrons is retained in soft Wilson lines, 5*^, with direction r^. 



22 



Since 5*^5'.,. = 1 these Wilson lines often cancel, but for many of the power suppressed terms 
listed in Table | the cancellation is not complete. This mechanism for generating a strong 
phase was first observed for D^tt^ 0, where a nonperturbative soft matrix element 

occurs through four-quark operators depending on n and v' (which are null and time-like 
vectors for the final state light and charmed mesons, respectively). 

For the B — > M1M2 decays with two energetic light mesons, a nonperturbative strong 
phase requires a soft matrix element depending on the Sn and Sn Wilson lines in SCETn. 
The simplest way to obtain the Wilson lines for the soft operators is to match SCETi onto 



SCETn P3. In SCETj one first uses the decoupling field redefinition on colhnear fields [|T6 

Yn^n, Yn^n, ^ YnAnY^ and An YnAfXn^ which generates the Wilson lines 

and factorizes usoft and colhnear fields. The fields of a given type are then grouped together 
by Fierz rearrangements. Matching the resulting operators or time-ordered products onto 
SCETn gives Yr — > Sr-, and we can read off which soft Wilson lines are present. Because 
of the properties of the subleading SCETi operators, we will not have an Sn and in the 
final SCETn operator unless we have a subleading SCETj Lagrangian with an n-collinear 
field and usoft fields, and one with n-collinear fields and usoft fields. We used this property 
to determine which entries are real or complex, and listed the results in the last column of 
Table |. The complex entries with multiple c!fq^^ 0] also have at least two hard-collinear 
gluons, and so generate contributions that start at as{^i)'^ when matched onto SCETn. 

To determine the perturbative order of the complex contributions, we must also classify 
which hard and jet coefficients give complex phases. In general any hard coefficient generated 
by matching at > 1 loop will give imaginary contributions, since these loops involve fields for 
both final state mesons, as pointed out for the general case in Ref. and for charm loops in 
Ref. [^. Since all leading order contributions in Table | have at least one ^^(/ij), the hard 
imaginary contributions for A^^^ are (9[as(/ii)as(/i/i)/7r]. At order K/mi, all annihilation 
contributions but Q^^ have at least one asijJ^i), and for these terms the hard complex 
contributions involve as{^ii)as{^ih) and thus are smaller than the nonperturbative terms 
proportional to Q;s(/ii)^. For Qf"* the amplitude is real at the leading perturbative order, 
as{fih), as demonstrated in section |T|, and so hard complex contributions start at 
In contrast for the amphtude A^^^^ a complex amplitude is generated at order as{Hi) A/ nib, 
which is only suppressed by A /nib compared to A'^^\ 

Finally, we should examine complex contributions from the jet scale. At leading order 
there is a unique jet function J 0. J also contributes to the heavy-to-light form factors and 
only knows about the n-collinear direction. Thus A^^^ does not get imaginary contributions 
at any order in the as{fii) expansion (which has been demonstrated explicitly to a^fii) [p8[]). 
At next-to-leading order in the power expansion, there is no known relation of the power 
suppressed jet functions with analogous jet functions in the form factors. However, the 
subleading jet functions also depend only on one colhnear direction, and do not carry in- 
formation about both final state mesons that could generate a physical strong phase. We 
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FIG. 3: Graphs which generate a strong phase in lowest order matching of SCETj operators onto 
SCETtt: a) has a Q^^\ two c9'^ , and one c9^^ and contributes to the annihilation amplitude 

at Oia^JLii)): and b) has a Q^^' , one C''}^ „, and one vCl^^ c and contributes to non-annihilation 
amplitudes at C'(as(/Xj)). Dashed quark lines are n\ or n2 collinear, and solid quark lines are soft. 



demonstrate this fact more explicitly by examining the calculation at (^(^^(yUj)), which is 
sufficient to see that the amplitudes are real up to the order where a nonperturbative phase 
first occurs. At this order the jet functions are generated by matching tree level SCETi 
diagrams onto SCETn. A typical example is 

^ (48) 



[x + ze) + ie 



where x is a momentum fraction that will be convolved with a collinear distribution function, 
and the k'^ will be convolved with a soft distribution function. These jet functions are real 
if and only if we can drop the ie factors. However, just as in section |T^, the ie terms 
can be dropped because the zero-bin subtractions |]18[ ensure that this does not change the 
convolution.'' Thus factorization gives real 0{as{\ii)) jet functions. 

This demonstrates that complex contributions in the power suppressed annihilation am- 
plitudes are suppressed. 



Im 



On general grounds one might have expected 0{A/mf,) suppressed strong phases, which we 
have demonstrated are absent in ^inn, though they do occur in A^^J^f^. 

We close this section by giving two examples of time-ordered products generating the 
nonperturbative strong phases discussed above. We consider a time-ordered product with 
three L^^^ insertions contributing to annihilation. When matching onto SCETn we integrate 
out the hard-coUinear modes, leading to an eight-quark operator. Figure ^ shows the order 
a'^ifii) contribution to this matching. The soft quark lines remain open as their contraction 



A equivalent physical argument for dropping the ie factors was given in Ref. js), where it was needed to 
prove that certain long-distance contributions are absent in color suppressed decays. 
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leads to an on-shell line which must be treated nonpertrubatively. The resulting SCETn 
operator has the generic form 

O^^ = J{n2 ■ p,ni ■ I, rii ■ r, n2 ■ g, ni ■ k) (50) 

where we use the shorthand subscript notation, (S'^j,gs)„..g = [6{ni-q — ni-V)Sl_qs]. We took 
the jet directions to be ni and n2, rather than n and n, to emphasize that the soft operator is 
sensitive to the relative directions of the jets. The functions Si shown in Table | are defined 
by the matrix element of this type of operator 

Si{ni-k,ni-r,n2-q,) = {0\{qsSn,)ni-r'r[^\si^qs)n2-q {qsSn2)n^-k'i^? {Si^K)\B{v)) , (51) 

where i runs over color, Dirac, and flavor structures. To count the factors of vr in these 
amplitudes, note that the hard-coUinear contractions give g'^, and that the matrix element 
of the resulting four-quark operator, {0\{q . . .q){q . . .by)\B), is suppressed by l/(47r)^ rel- 
ative to {0\{q . . . bv)\B) . (The four-quark operator has an extra loop with no extra cou- 
plings.) This demonstrates that nonperturbative complex contributions first occur at order 
[as(/ii)^/7r](A/m;,), i.e., suppressed by [as(yUj)/7r](A/m;,) compared to the leading amplitudes. 
The phases arising from the type of matrix element shown in Eq. ( pl| ) play a crucial role 
in explaining the observed strong phases which arise in color suppressed decays 0. Their 
resulting operators predict the equality of amplitudes and strong phases between decays 
involving D and D* mesons and have been confirmed in the data This type of dia- 
grams also have long-distance contributions of the same order, which arise from time-ordered 
products in SCETn and can also be complex. To see this note that the hard-coUinear quark 
propagator in Fig. |^a could also be on-shell (i.e., have (9(A^) virtuality), in which case it 
would remain open until the matrix element is taken at the low scale. By opening that line 
we see that this contribution corresponds to the time-ordered product of a four-quark op- 
erator and a six-quark operator, both of which are generated when matching onto SCETn. 
A long-distance part is the same order in ^^(/ij) and does not change our conclusions about 
these terms. In Fig. ^ we show a non-annihilation contribution to A^g^^ which is of order 
as{f^i)A/mh. This term is generated by the time-ordered product of Q^^\ an insertion of the 
rii-coUinear jC^^g , and an operator with n2-collinear quarks and usoft gluons. 



Cf^ = (LW) Ylilj^, ip^^^Y^, £ {W^i,) . (52) 



VI. APPLICATIONS AND CONCLUSION 
A. Phenomenological Implications 

To understand the implications of the experimental data, it is crucial to know which 
contributions to the B M1M2 amplitudes can be complex. The best sensitivity to non- 
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SM physics is via interference phenomena, where new interactions enter hnearly (instead of 
quadratically) , such as CP- violating observables. The sensitivity to such effects depends 
on how well we understand the dominant and subdominant SM amplitudes, including their 
strong phases. The existence of strong phases in B decays is experimentally well established 
(e.g., the B Dn and S — > vrvr rates, the CP asymmetry Ax+n-, the transversity analysis 
in B J/ipK*, etc.). 

One example of how strong phase information can be useful is the method for determining 
7 from B —* TTTT proposed in Ref. The method uses isospin, the factorization prediction 
that Im(C/T) ~ 0(a,(m;,) , A/nih) , and does not require data on the poorly measured direct 
CP asymmetry Ct^Oj^q.^ The phases in A^^^ at Q;s(mf,)as(/ij) are calculable and partially 
known p|, ^ . The current B ^ mr data is in mild conflict (at the ~ 2a level) with the SM 



CKM fit More precise measurements are needed to understand how well the theoretical 
expectations are satisfied, and to decipher whether there might be a hint for new physics. 
Obviously further information about power corrections in Im(C/T) could help to clarify the 
situation. 

In all factorization-based approaches to charmless B decays, several parameters are fit 
from the data or are allowed to vary in certain ranges. The choice and ranges of these 
parameters should be determined by the power counting. This motivated keeping the charm 



penguin amphtudes, Acc as free parameters in SCET P, as was done earlier in Ref. []T2 
In the BBNS approach these are argued to be factorizable 0. A fit to the data using this 
parameterization found large power suppressed effects |^ including annihilation amplitudes, 
which might be interpreted as a breakdown of the A/mb expansion. In QCD sum rules, the 
annihilation amplitude was found to be of the expected magnitude and to have a sizable 
strong phase [Q, but a distinction between the terms we identify as real local annihilation 
and complex time-ordered product annihilation was not made. 

Channels like B Kn and B — ^ KK are sensitive to new physics, but by the same to- 
ken are dominated by penguin amplitudes, which can have charm penguin, annihilation, and 
other standard model contributions. Since there are possible large nonperturbative c-loop 
contributions in Acc that have the same SU (3) flavor transformation properties as annihi- 
lation terms, they cannot be easily distinguished by simple fits to the data. However, in a 
systematic analysis based on SCET these correspond to different operators' matrix elements, 
so it is possible to disentangle the various contributions and determine their expected size. 
The factorization theorems for annihilation amplitudes derived here only involve distribu- 
tions that already occurred at leading order. This means that we can compare the size of 
annihilation amplitudes to experimental data without further ambiguities from additional 
hadronic parameters. We take up this comparison in section |VlIj| below. 



^ Here C and T are isospin amplitudes defined in the i-convention, where At is ehminated from the amph- 
tudes in favor of Ac and Au. 
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As an explicit example of how to assemble our results in sections |Tl| and |I^, we derive 
the local annihilation amplitude for K~tt~^. From Table ^ we can read off the result 

for this channel, H{x,y) = —al{x,y) — ag(x,?/), and from Table H^i = —a^{x,y) + 
1/2 a^{x,y) and H^2 = 0-2 (3^,2/) — 1/2 aQ(a;, y). With the lowest order matching results in 
Eqs. (El) and (P^) we can set ag = and a^u = (^ic, which inserted into Eqs. and 



gives 
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Thus, both the leading order annihilation amplitude A)^/^^^, and the chirally enhanced anni- 
hilation amplitude v4^^j^„ are determined by the /3's defined in Eqs. ( p6|) and (^i]). Other Ktt 
channels have similar expressions with different Clebsch-Gordan coefficients. To the local 
annihilation contributions we must add the hard-coUinear annihilation terms computed in 
Ref. |]19|, ^hard-coiiin' siucc they are the same order in as and l/rrib as the A^^^^ terms. 
To see exphcitly what the /3's involve we insert the 0{as) values of a3u{x,y), a^{x,y), and 
(x, y) to give 
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+ d,{f,_)Mi){{y~T+{y~T)-d3if^±)<P'MM^) 



Here results for the convolutions denoted by brackets (■ ■ ■ ) can be found in Eqs. (p9D, 
(0), (|46|), and (|47|) in the minimal subtraction scheme. Results for other channels can 



be assembled in a similar fashion. Corrections to A^ann + ^hird-coUin suppressed by 
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0[a'^{fii) / {nas{mh))], while we caution that additional as{fih)-^/mb terms without a /i,r or 
fix will be present in the last two lines of Eq. (p^. In the next subsection we derive results 
for all of these channels using a simple model for the distribution functions, and study 
numerically the size of the annihilation amplitudes. 

Annihilation contributions have been claimed to play important roles in several observ- 
ables H, |1^, [Tl|, in particular in generating large strong phases in B ^ Kn de- 



cays [0, |]. The B ^ nil and Kn data indicate that the latter decays are dominated by 
penguin amplitudes, and the pattern of rates and CP asymmetries is not in good agreement 
with some predictions. In particular, it is not easy in the BBNS analysis to accommo- 
date the measured CP asymmetry, Ak+-k- = —0.108 ± 0.017 Esl, except in the S3 and 



S4 models of Ref. In these models the annihilation contributions are included by us- 



ing asymptotic distributions, and divergent integrals are parameterized as J^dx/x — > 



and J^dxlnx/x -Xl/2, with Xa = {I + ^Ae**^-^) ln(mB/500 MeV). Model S3 pos- 
tulates qa = ^, '^A = —45° for all final states, while in the S4 scenario Qa = ^ and 
= -55°, -20°, -70° for the PP, PV, VP channels, respectively. Thus 

53 : = 4.0 - 1.7i, 54 : = {3.7 - 1.9 z , 4.6 - 0.8 z , 3.2 - 2.2 z} . (55) 



In addition, as(/i) and the Wilson coefficients are evaluated at the /Xj intermediate scale |]TT|. 

Our result for the factorization of annihilation contributions derived in Sec. |ITT| constrains 
models of annihilation. Equation ( |23| ) gives a well defined and real amplitude at leading 
order, which depends on twist-2 distributions, 0m- It does not involve model parameters qa 
and (fiA- For A^H^^ using Eq. ( P^ and the asymptotic form of the meson distributions, we 
find a correspondence 

'•^-." = 1+ A- 5^ = '"(-)■ (56) 

Jo Q[x^)0 

Clearly, Xa is real. The asymptotic distributions ~ 6x(l — x) are more accurate for large 
scales, and at the matching scale where /i+ ~ m^, Xa is not enhanced by a large logarithm. 
This matching scale /i+ should not be decreased below since /i+ ~ is already the 
correct scale for collinear modes with ~ rrib- We estimate \Xa\ ^ 1. Thus, the modeling 
of annihilation contributions with complex Xa in the BBNS approach (including the phe- 
nomenologically favored S3 and S4 scenarios) are in confiict with the heavy quark limit, and 
should be constrained to give smaller real X^'s. 

In the KLS treatment of annihilation, complex amplitudes are generated from dy- 
namics at the intermediate scale from the ie in propagators. The MS-factorization used in 
the derivation of our annihilation amplitudes demonstrates that including the ie term in 
collinear factorization would induce a double counting. Thus we expect such contributions 
to physical strong phases to be realized by operators with soft exchange that occur at higher 
order in A/rrib and therefore to be small. 

Annihilation contributions were also argued to play an important role in explaining 
the large transverse polarization fraction in i? ^ (f)K* |^. It was shown that fac- 
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torization implies Rt = 0{l/ml), where Rt denotes the transverse polarization frac- 



tion |]31]]. Subsequently, it was shown using SCET that Rt is power suppressed unless a 
long-distance charm penguin amplitude Acs spoils this result Experimentally, one 

finds Rt{B —>■ (f)K*) ^ 0.5 HH], while Rt{B —>■ pp) is at the few percent level. It has been 



argued that the large Rt{B (pK*) may provide a hint of new physics in the b sss 



channel. In Ref. |^T| it was suggested that standard model annihilation contributions may 
account for the observed large value of Rt{B — > (pK*). Our analysis in Sec. |V| agrees 
with pri in that annihilation contributions to the transverse polarization amplitude at first 



order in are suppressed by not one, but two powers of kjm^,. However, we do not find 
a numerical enhancement of these terms (which in ||3l| is partly due to the large sensitivity 
of the (2X^ — 3)(1 — Xa) function to qa in the BBNS parameterization). The operators in 
Eq. (|36D give rise to transverse polarization, but since MS-factorization renders the naively 
divergent convolutions finite, these power suppressed amplitudes do not receive sizable en- 
hancements. Although we have not derived explicit results for the B —>■ (j)K* annihilation 
amplitudes (since is an isosinglet), our results make it unlikely that local annihilation can 
explain the Rt{B <f)K*) data. We have not explored whether the time-ordered products 
at 0{a'^{fii)A/'mb) could give rise to transverse polarization, and it would be interesting to 
do so. 

B. Annihilation amplitudes with simple models for (f)^{x) and </>^(ic) 

In this section we derive numerical results for the local annihilation amplitudes in various 
channels using a simple model for the distributions. It is convenient to write the AS" = 
local annihilation amplitude as 

ALanniB ^ M^M^) = -^iM^Ml i^X^^)K{B ^ M.hh) + X^^ K{B M.M,) 

+ (Ai') + Ai'^)) h^,{B ^ M1M2) + ^ h^2{B ^ MiMa)] \ . (57) 
I nib rnb -I J 

For AS* = 1 decays we replace \u}c — > \u}c- The coefficients h^, he-, h^i, and h^2 are equal 
to linear combinations of /3jc, /^xi? /3x2' /^xS' ^"^^ /^x6 with Clebsch-Gordan coefficients 
determined from Tables y, [m] , |V|, |V|. The combinations are simply determined by the 
replacements 

K = {H{x,y) with ~af\x,y) - df\y,x) - /3f^^^^) , 

K = {H{x,y) with ~af%x,y) - ^^ df\y,x) - *^^) , 
^xi = {Hxiix,y) with alr^{x,y) Pxhx^) > 

^x2 = {Hx2{x, y) with a|6(a;, y) (3x2,x6) ■ (58) 
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For the coefficients ci3u,3c,4u,4c,7,8 and the af^s, the C(a^Ci^2) matching corrections could be 
comparable numerically with the OlagC^^io) corrections considered here. This should be 
kept in mind when examining numbers quoted below for the corresponding /5's. 

Results for the coefficients Piu, Pic, and P^i, can be found in Eqs. and (|^). To derive 
numerical results we need to model the meson distribution functions. We take the Cj from 
Eq. use 

asifih) = 0.22 , fi^{fih) = 2.3GeV, fiKif^h) = 2.7 GeV , 

/K = 0.16GeV, /^ = 0.13GeV, /B = 0.22GeV, (59) 

where /ih = '"^6 = 4.7 GeV, fs comes from a recent lattice determination [Q. For the 0's 
we take simple models with parameters a*^ and af^^ which we consider specified at the high 
scale fih, 

^^\x) = 6x{l -x)[l + af (6x - 3) + 6af (1 - 5x + 5x^)] , 



6x(l — x) [l + a{^p(6x — 3) + 6a^p(l — 5x + 5x^)] . (60) 



Based on recent lattice data for moments of the vr and K distributions |57| we take 03 



0.2 ± 0.2, where the lattice error was doubled to give some estimate for higher moments. 
For the vr we set = a^^^ = 0, while for the K we use |^ af = 0.05 ± 0.02. We also take 
wz-K^K = — 3 ± 1, Ogp^ = ± 0.4 and a{^p = 0.0 ± 0.2. Note that the range for our parameters 



is similar to those used in the BBNS models |Ty, |TT| and light-cone sum rules . Since the 
uncertainties in the model parameters are large and not significantly affected by variation of 
the /x± scales we keep these fixed at m^, where the logs in the di{fi±) terms drop out and the 
constant under the logs are neglected. A scan over models with parameters in these limits 
gives predictions for the annihilation coefficients. For the B —>■ Kn channels we find 

= 1.8 ± 1.2, = = -0.15 ±0.10, = 0.14 ± 0.09 , 

Pl^^ = 0.09 ± 0.33, pl^2 = -0-29 ± 0.09, = -0.012 ± 0.002, = 0.002 ± 0.01 , 

p^^ = 0.0 ±6.5, Pl^ = 0.0 ±5.8, = 0.0 ± 0.094, = 0.0 ± 0.11 . 

(61) 

Using these numbers we can compare the size of the local annihilation amplitudes to the 
B — >■ K~7i^ data, 

R^{K-n+) = I^Snn(-^"^+) + = Q.ll ± 0.09 , 

R^{K%-) = + ^ Q^^2 ± 0.09 . (62) 

\^Expt.Penguin[-i^'^)\ 

For the numerator we did a Gaussian scan using the values from Eq. (0), and determined 
the error by the standard deviation. For the denominator we used the experimental penguin 
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amplitude determined by a fit to the B Kn data in Ref. B. Numerical results for 



annihilation amplitudes with three-body distribution functions were considered in Ref. p9 
Although they are similar in size to ^^^i„„ they cause only a ~ 10% change in the value of 
Ryi{K~7r'^) in Eq. (^). The values of Ra indicate that a fairly small portion of the measured 
penguin amplitude is from annihilation. We do not quote values for the ratio A^Lann/^Lanni 
since each of the numerator and denominator can vanish and the parametric uncertainties 
are very large. For typical values of the parameters in the Ktt channels we find that the 



^Lann Comparable or even larger than A^H^^ in agreement with Ref. |T^. The size of 



the annihilation amplitudes in Eq. (|62D are consistent with our expectation for these power 
corrections. For B KK we find 

(31'' = -9.6 ± 6.2, = 1.7 ± 1.1, = Pi/ = 0.63 ± 0.37, 

Pij< = p^K ^ _o.i4 ± 0.09, P^/ = -0.03 ± 0.02, p^/ = 0.13 ± 0.08, 

p^Ji = pg^ = 0.63 ± 0.37, p^f = 0.0 ± 6.5, p^/ = 0.0 ± 5.5 

/3j^ = 0.0 ± 0.095, = 0.0 ±0.11. (63) 

Using these results to determine the Xi^"^ annihilation contributions to -B ^ KK and com- 
paring this to the experimental penguin amplitude from Ref. gives 

RAiK~K') = \^^Lann{^^K^) + AfanniK' K^)\ ^ ^ ^ ^^^^ 

\AExpt.Penguin[KK )\ 

This is similar in size to the ratios Ra{K~7;~^), Ra{K^tt~) and so also consistent with a 
power correction. 



C. Conclusions 



In summary, we exhibited how a new factorization in SCET renders the annihilation and 
"chirally enhanced" annihilation contributions finite in charmless nonleptonic B M1M2 
decays to non-isosinglet mesons. We constructed a complete basis of SCETn operators for 
local annihilation contributions as well as factorization theorems valid to all orders in as- 
By matching the full QCD diagrams onto SCETn operators we showed that their matrix 
elements are real at leading order in A/rrif, and asifrif,). The lowest order annihilation contri- 
butions depend on fs and a modified type of twist-2 distributions cj)^'^'^ with dependence on 
rapidity cutoffs. Chirally enhanced local annihilation contributions depend in addition on 
modified distributions 0^^'^. The annihilation contributions can only have an unsuppressed 
complex part at 0{A/mh) if perturbation theory at the intermediate scale, ^/Amb, breaks 
down. 

In the previous literature models for the power suppressed annihilation corrections were 
often found to give enhanced contributions with large strong phases, and such assumptions 
have been important in some fits to the data. Considering all power suppressed amplitudes 
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not involving charm loops, we proved that complex annihilation contributions only occur 
suppressed by as(v^Amfe) Agco/'^b compared to the leading amplitudes. From our factor- 
ization theorem we found that annihilation contributes (11 ± 9)% of the penguin amplitude 
in 5° ^ (12 ±9)% in 5" i?°7r-, and (15 ±11)% in K'K^. We anticipate 

that our results will guide future fits to the vast amount of data on charmless B decays, and 
yield a better understanding of what this data means. 
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APPENDIX A: ZERO-BIN SUBTRACTIONS FOR A TWO-DIMENSIONAL 
DISTRIBUTION 



In this appendix we derive a result for the action of the zero-bin subtractions on the 
integrand obtained from the chirally enhanced annihilation computation, shown in Eq. (^7|). 
Since the result involves a correlation in the x and y integrals it cannot be read off from the 



results in Ref. [0. It is convenient to write the momentum fraction factor coming from the 
offshell 6-quark propagator as (1 — xy) = {x + y — xy). Including the rapidity convergence 
factors |T^, the integral we need is 



where O^^ = 6{x)6{l — x). To determine the subtraction terms we must look at the singular 
behavior as we scale towards the x = 1 and y = bins, which we do by taking x ^ rj and 
y ^ 1]. In this limit the gluon and 6-quark in Fig. |^ become soft, and this region would be 
double counted without the zero-bin conditions. First consider the denominator, 

1 I xy , . 

+ 7— + • • • • (A2) 



X + y — xy {x + y) {x + y) 

In the first term the x and y dependence does not decouple, so we must consider them 
simultaneously. All terms beyond the first one produce finite integrals and are dropped in 
the minimal subtraction scheme. For the numerator in Eq. ( [Al| ) we use 4>^^{0) = (/)^^(1) = 
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and expand 



#;,(0) - ^C,(0) a;0'(l) + #;,(0) - r (1) + . . . 



VP 



(o)x0'(i) + ... . 



n=l 



n! 



(A3) 



In the first term on the last hne we have identified all terms which remain singular when 
multiplied by l/[xy'^{x + y)]. This term is equal to ?/0pp(O)0(x). Taken together with the 
expansion of Qx^y we therefore find that the required minimal subtraction is 



(A4) 



{x + y)xy'^ 

Following Ref. [1^ we use this to convert Eq. ( |A1|) into an integral that includes the x = 1 
and y = regions, 



dx I dy 



1 

dx 

sM2 



{x + y — xy)xy'^ 



{x + y)xy'^ 



dy \ ^^(1 - y\y - ^ 



{x + y)xy 



2e 



(A5) 



dx ■ 



x] 



X 

kM2l 



dx ■ 



X) 



X 



dy 



dy 



_ 

{x + y- xy)y'^ (x + y)y 



dx 



C^'(O) / dx 



{x + 

ln(2-x) 



Here in simplifying the term carrying the y ^ oo limit, we noted that the integral is finite, 
and so it does not induce fi± dependence in our subtraction scheme. This result for / 
was used in Eq. (0). For the asymptotic pion wave functions, (f)'^{x) = 6x{l — x) and 
(j)^p{y) = 6y{l - y), we obtain J = 36 + Ovr^ - 144 In 2 = -4.60. Note that the steps used 
here to derive the subtraction also give the correct result for cases where the x and y integrals 
factorize, such as an integrand (f){x)(f){y)/ {x'^y'^). 
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